Describing a First order System Using Matrix Notation 

1. Description of the Equation 

A general 2x2 linear system is given by: 

X = ax + by 
y = cx-\- dy 

The terms have been arranged in a suggestive manner. We can express this 
system using matrices and vectors: 

X \ _ ( a h \ f X 
y ) ~ \ c d )\y 

We can present this in the following even more compact form. 

i u h \ IX 
Let A = I , and write u for the colimm vector 

\ c d J \y 

We have u(t) = ( ^^^j 1 and the system is simply u = Au. 

Example 1. Our favorite system, governing the rabbit populations in farm- 
ers Jones' and McGregor's fields, was 

X = 0.3x + O.ly 
y = 0.2x + 0.4y, 

which has matrix form 

x \ / 0.3 0.1 \ / X \ . „ , „ / 0.3 0.1 

y ) = [o.2 OA)[y ) " = "^^^'^ ^ = [ 0.2 OA 



2. Description of the Solution 

To describe the solution, we will use the column vector u{t) 



x{t) 

y{t) 



Example 2. Earlier we used the method of elimination to solve the system 
in example 1. We found x{t) = cie^-^* + Cie"'^*, y{t) = 2cie°-^* - C2e°-^K 
Rewriting this in vector form we have 



u(0 
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We can rewrite this as 

which is a clearer way of presenting it. Let 

Mt)=e''[l^ andu2(0=.-^^(_J). 

The colimrn vectors ui(t) and U2{t) are both solutions. Since they both 
involve only one form of exponential, they are sometimes known as basic 
independent solutions, or normal modes. The general solution is a linear 
combination of them. We will learn much more about normal modes in the 
sessions on matrix methods and the phase portrait. 

Remark. As with linear second order ODE's in unit 2, the general solution 
to to a 2 linear system should always consist of linear combinations of two 
truly different solutions. It is not necessary , but usually our techniques 
will make these two solutions the normal modes. 

3. Geometry of the Solutions 

Suppose you want to plot a solution u(f). As time increases, it traces a 
curve in the xy-plane. 

Example. The solution ui{t) traces a ray that passes through (1,2) at i = 
and move off towards infinity in a straight line, with exponential speed. 
There is another ray through (1, —1), corresponding to the solution U2{t). 
(This is tricky: the exponential in the formula for ui might make you think 
the trajectory is curved. However, if you look carefully at the formula you 
will see ui (t) is always a multiple of the vector (1,2)^.) 

The applet Linear phase portrait: matrix entries will allow us to visualise 
this nicely, and get a feel for other sorts of trajectories. Later in this session 
you will look at this applet. 
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